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Abstract 

We search for lepton-flavor-violating r — )• iV^ decays, where i is an electron or muon and 
is one of the vector mesons p^, cj), w, K*^ and K*^. We use 854 fb~^ of data collected with the 
Belle detector at the KEKB asymmetric-energy e'^e~ collider. No evidence for a signal is found 
in any decay mode, and we obtain 90% confidence level upper limits on the individual branching 
fractions in the range (1.2 — 8.4) x 10~^. 

PACS numbers: 11.30.Fs; 13.35.Dx; 14.60.Fg 



3 



INTRODUCTION 



Lepton flavor violation (LFV) in charged lepton decays is forbidden in the Standard 
Model (SM) and remains highly suppressed even if the SM is modified to include neutrino 
mixing. However, extensions of the SM, such as supersymmetry, leptoquark and many other 
models [l|-0] predict LFV with branching fractions as high as 10~^, which could already be 
accessible in current 5-factory experiments. We search for r~ — ?■ £~V^ decays [f], where i is 
an electron or muon and is one of the vector mesons p°, 0, u, K*^ or K*^. These results 
are based on the entire 854 fb~^ data sample collected at center-of-mass (CM) energies near 
the T(4S'), near the T{5S), and off resonance with the Belle detector [sj at the KEKB 
asymmetric-energy e'^e~ collider j^. Previously, we obtained 90% confidence level (CL.) 
upper limits on branching fractions of these decays using 543 fb~^ of data; the results 
were in the range (5.8—18) x 10~* [13] • The results described here use additional data 
and an improved event selection procedure, which is optimized mode-by-mode. The BaBar 
collaboration has also published 90% CL. upper limits in the range (2.6—19) x 10~^ using 
451 fb-i of dataO for all r i decays except for r — )■ £ a; for which 384 fb ^ of 
data were used 12 1. 



The Belle detector is a large-solid-angle magnetic spectrometer that consists of a silicon 
vertex detector (SVD), a 50-layer central drift chamber (CDC), an array of aerogel threshold 
Cherenkov counters (ACC), a barrel-like arrangement of time-of- flight scintillation counters 
(TOP), and an electromagnetic calorimeter comprised of CsI(Tl) crystals (ECL), all located 
inside a superconducting solenoid coil that provides a 1.5 T magnetic field. An iron flux- 
return located outside of the coil is instrumented to detect mesons and to identify muons 
(KLM). The detector is described in detail elsewhere [8|. 

Particle identification is very important for this measurement. We use particle identi- 
fication likelihood variables based on the ratio of the energy deposited in the ECL to the 
momentum measured in the SVD and CDC, shower shape in the ECL, the particle's range 
in the KLM, hit information from the ACC, dE/dx measured in the CDC, and the particle's 
time of flight. To distinguish hadron species, we use likelihood ratios, V{i/ j) = Ci/ (Ci + Cj), 
where Ci (Cj) is the likelihood of the observed detector response for a track with flavor i 
(j). Por lepton identification, we form likelihood ratios V{e) [l^ and V{n) [13] using the 
responses of the appropriate subdetectors. 

We use Monte Carlo (MC) samples to estimate the signal efficiency and optimize the 
event selection. Signal and background events from generic t~^t~ decays are generated by 
KKMC/TAUOLA [l^. Por signal, we generate e"'"e~ — )■ t~^t~ events, in which one r is 
forced to decay into a lepton and a vector meson using a two-body phase space model, 
while the other r decays generically. Background events BB, continuum e'^e~ — )■ qq 
u,d,s,c), Bhabha, and two-photon processes are generated by EvtGen [l^, BHLUMI 
and AAPH [3], respectively. In what follows, all kinematic variables are calculated in the 
laboratory frame unless otherwise specified. In particular, variables calculated in the e'^e~ 
CM frame are indicated by the superscript "CM" . 



[t] Throughout this paper, charge-conjugate modes are imphed unless stated otherwise. 
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TABLE I: Allowed ranges for the mass of V candidates and for the cosine of the angle between 
the missing particle and the tag-side track in the CM system (cos 





Invariant mass (GeV/c^) 


cosC^„^issforr^^yO (eV^) 




0.587 < M^^ < 0.962 
1.009 < Mkk < 1-031 
0.757 < M^^^ < 0.808 
0.842 < Mxn < 0.956 


[0.0, 0.85] ([0.0, 0.96]) 
[0.0, 0.88] ([0.0, 0.97]) 
[0.0, 0.88] ([0.0, 0.97]) 
[0.0, 0.87] ([0.0, 0.96]) 



EVENT SELECTION 

We search for e~^e~ — ?■ r^r~ events in which one r (the signal r) decays into a lepton 
and a vector meson [V'^ = p^, 0, u, K*^ and K*^) and the other r (the tag r) results in 
a single charged track, any number of additional photons, and neutrinos. We reconstruct 
p° from TT+TT", (j) from K~^K~ , u from tt^tt'tt^ , K*^ from K^'k~ and K*^ from K~tt^ . To 
improve sensitivity compared to our previous work [l^, the event selection is optimized 
mode-by-mode, since the backgrounds are mode dependent. 

For each candidate event we calculate the iV^ invariant mass (M^v") and the difference 
of the IV^ energy from the beam energy in the CM frame (Ai?). In the two-dimensional 
distribution of M^yo versus Ai?, signal events should have M^yo close to the r-lepton mass 
(m,-) and Ai? close to zero. 

We require that events contain four charged tracks and any number of photons within the 
fiducial volume defined by —0.866 < cos^ < 0.956, where 9 is the polar angle relative to the 
direction opposite to that of the incident beam in the laboratory frame. The transverse 
momentum {pt) of each charged track and the energy of each photon {E^) are required to 
satisfy Pi > 0.1 GeV/c and > 0.1 GeV, respectively. For each charged track, the distance 
of the closest point with respect to the interaction point is required to be less than 0.5 cm 
in the transverse direction and less than 3.0 cm in the longitudinal direction. 

Using the plane perpendicular to the CM thrust axis [l9|, which is calculated from the 
observed tracks and photon candidates, we separate the particles in an event into two hemi- 
spheres. These are referred to as the signal and tag sides. The tag side is required to 
contain one charged track while the signal side is required to contain three charged tracks. 
We require one charged track on the signal side to be identified as a lepton. Electron candi- 
dates are required to have P(e) > 0.9 and momenta p > 0.6 GeV/c while muon candidates 
should have V{fi) > 0.95 and p > 1.0 GeV/c. In order to take into account the emission of 
bremsstrahlung photons from electrons, the momentum of the electron track is reconstructed 
by adding to it the momentum of every photon within 0.05 radians of the track direction. 
With these requirements, the electron (muon) identification efficiency is 91% (85%) while 
the probability to misidentify a pion as an electron (a muon) is below 0.5% (2%). 

Charged kaons are identified by the condition V{K/tt) > 0.8 for the r — )■ iK*^ and iK*'^ 
modes (> 0.6 for the £(f) modes). Charged pions must satisfy the requirement V{K/ii) < 0.6 
for the r — ?■ CK*^ and iK*^ modes. We do not require charged pion identification for the 
T ^ £uj modes. The kaon (pion) identification efficiency is 80% (88%) while the probability 
to misidentify a pion (kaon) as a kaon (a pion) is below 10% (12%). In order to reduce the 
fake vector meson background from photon conversions (i.e., 7 — )■ e~^e~), we require that 
the two charged tracks from the meson candidate have V{e) < 0.1. Furthermore, we require 
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FIG. 1: Vector meson invariant mass distributions for the r — t- fiV modes after the basehne 
selection (before vector meson mass requirements) in the region 1.5 GeV/c^ < M^yo < 2.0 GeV/c^ 
and -0.5 GeV < AE < 0.5 GeV: (a), (b), (c), (d) and (e) show the r /i(/>, fiu, fiK*° and 

fj,K*^ modes, respectively. The signal MC (r — )■ ij,V^) distributions are normalized arbitrarily while 
the background MC distributions are normalized to the data luminosity. Selected mass regions are 
indicated by arrows. 



< 0.1 to suppress e~^e~ — )■ e~^e~ fi^fi^ two-photon background. For the r — )■ modes, 
we select 7r°(— 77) candidates on the signal side by requiring > 0.1 GeV and Pt,o > 0.4 
GeV/c, which suppresses incorrectly reconstructed ir^ candidates from initial-state radiation 
(ISR) and beam background. The vr" mass window is 0.12 GeV/c^ < M^^ < 0.15 GeV/c^. 

To ensure that the missing particles are neutrinos rather than photons or charged par- 
ticles that pass outside the detector acceptance, we impose requirements on the missing 
momentum jViss; which is calculated by subtracting the vector sum of the momenta of all 
tracks and photons from the sum of the e"*" and e~ beam momenta. We require that |Pmissl; 
the magnitude of the transverse component of Pmiss, be greater than 0.5 GeV/c (0.7 GeV/c) 
for the r — >■ fiV^ {eV'^ except the ep°) modes, and that Pmiss point into the fiducial volume 
of the detector. For the r — ep° mode only, we apply the more restrictive selection re- 
quirement |p|nissl > 1.5 GeV/c. Furthermore, we reject events if the direction of the missing 
momentum traverses the gap between the barrel and the endcap since an undetected photon 
may be incorrectly attributed to missing particles. 
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FIG. 2: (a) V{K/Tt) for muon tracks, (b) V^p/tt) and (c) V{p/K) for hadronic tags, for /ip^ candi- 
date events satisfying 1.5 GeV/c^ < M^yo < 2.0 GeV/c^ and -0.5 GeV < AE < 0.5 GeV. Signal 
MC (r — )■ fip^) distributions are normalized arbitrarily while the background MC distributions are 
normalized to the data luminosity. Selected regions are indicated by arrows. 



To suppress BB and qq background, we require that the number of photons on the tag 
side, n^^'^, satisfy n^^*^ < 2 (< 1) for hadronic (leptonic) tags. For all modes except 
71-0 _j, g^jgQ rgq^irg at most one photon on the signal side. 

We select vector mesons whose invariant mass satisfies the requirements shown in Table |T] 
and Fig. [1] MC simulation predicts that for the r — )■ fiV'^ modes the dominant background 
comes from continuum and generic r"'"r~events, whereas for the r — )■ eV^ modes it originates 
from inelastic l^°-photoproduction {e~^e~ — )■ e~^e~V^) and two-photon processes. 

Since neutrinos are emitted on the tag side only, the direction of jViss should lie within 
the tag side of the event. The cosine of the opening angle between Pmiss and the charged 
track on the tag side in the CM system, cos 6'^j^g_j^jgg, is required to be greater than zero. 
If the track on the tag side is a hadron, then for the r — )■ ^V^ modes we also require 
cos ^J?j^g_n^jgg < (0.85 — 0.88). This cos^^^_j^jgg requirement reduces continuum background 
with missing energy due to neutrons or K^s since the masses of the neutron and Kl are 
higher than that of the neutrino. We also require that cos 6'^ag_j„igg < (0.96 — 0.97) for the 
r — )• eV^ modes. This requirement can reduce Bhabha, inelastic V°-photoproduction, and 
two-photon background, since radiated photons from the tag-side track result in missing 
momentum if they overlap with the ECL clusters associated with the tag-side track (20| . 
The cos 6'-^^'^ requirements for all modes are given in Table HI 

For the r — )■ fiV^ modes, a muon can originate from a kaon decaying in the CDC [K^ — )■ 
/i^z/^). Therefore, we apply a kaon veto, V{K/7t) < 0.6, for muon candidate tracks if the tag 
side track is a hadron (see Fig. [2] (a)). Another important continuum background predicted 
by MC is from dibaryon production with a proton on the tag side. Therefore, we apply a 
proton veto, V{p/n) < 0.6 and V{p/K) < 0.6, as shown in Figs. [2] (b) and (c). To reject qq 
continuum background, we require that the magnitude of thrust (T) be greater than 0.90. 
The reconstructed mass of the charged track and photons on the tag side, mtag, is required 
to be less than the r lepton mass (1.777 GeV/c^). To reduce e'^e~ background for 

the r -7- fip^ mode only, we also require that the momentum of the track on the tag side in 
the CM system be less than 4.0 GeV/c if the track is a muon candidate. 

For the r — )■ eV'^ modes, photon conversions can result in large backgrounds when an 
e~^e~ pair is incorrectly reconstructed as a fake vector meson. Fake vector mesons can also 
originate from one of the electrons in a photon conversion and a hadron. As seen from 
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FIG. 3: Scatter-plot of the momentum of kaon candidates from A'**^ candidates, Px+i versus 
the invariant mass of electron and kaon candidates from r — >■ eK*^{^ K~^ir~) reconstructed by 
assigning the electron mass to the kaon tracks, Mee{e~ K~^), for the eK*^ mode: (a), (b) and (c) 
show the distributions for data, generic rr MC and rr(-;> eK*'^) MC, respectively. Selected regions 
are indicated by red lines and arrows. 



generic t~^t~ MC, such a background comes from — h~Tx^{-^ 1'1)^t decays. When a 7 
from the vr^-decay converts and one of the conversion electrons traverses the ECL gap, then 
it may incorrectly be classified as a hadron because the electron veto does not work. In the 
momentum range between 0.5 GeV/c and 1.0 GeV/c such an electron can be misidentified as 
a kaon because the dE/dx distributions of kaons and electrons overlap. In other momentum 
ranges the electron can be misidentified as a pion. If we assign an electron mass to such a 
fake hadron, then the invariant mass Mee{e.~h^) will be small, as expected for conversions. 
Therefore we require Mee{e~ h'^){Mee{h~ h~^)) > 0.2 GeV/c^ for e~h~^ {h~h^) candidates for 
the ep*^, eK*'^ and eK*^ modes. Background from such e~K^ events in the eK*'^ mode is 
shown in Fig. [31 For the t ep^ mode only we also require that the magnitude of the 
thrust be in the range 0.90 < T < 0.96 to reduce two-photon and inelastic y°-production 
backgrounds. 



SIGNAL AND BACKGROUND ESTIMATION 

For all modes, the M^yo and AE resolutions are obtained from fits to the signal MC 
distributions, using an asymmetric Gaussian function that takes into account initial-state 
radiation. These Gaussians have widths as shown in Table [TTl 

To evaluate branching fractions, we use elliptical signal regions with major and minor 
axes equal to 3a in the AE and M^yo distributions, where cr = [a^^^^ + a^°'")/2 from Table [TTl 
The center and rotation angle are determined by scanning to maximize the efficiency of the 
signal divided by the area. We blind the data in the signal region until all selection criteria 
are finalized so as not to bias our choice of selection criteria. 

Figures [Hand [SI show scatter-plots for data events and signal MC samples distributed over 
±20(7 in the Miyo—AE plane for the r — ?• fiV'^ and eV'^ modes, respectively. For the r — >■ eV^ 
modes, the dominant background comes from two-photon processes, while the fraction of 
qq and generic t~^t~ events is small due to the low electron fake rate. For the r — )■ np^ 
mode, the dominant background comes from e~^e~ — )■ qq events. A smaller background of 
generic t~^t~ events in the region AE < GeV and M^yo < m-r is due to combinations 
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TABLE II: Summary of M^yo and AE resolutions {(Tm^'^ (MeV/c^) and cj^'f (MeV)). Here 
^high (^low) jneans the standard deviation on the higher (lower) side of the peak. 



Mode 






high 
"AE 


^low 
^AE 


T - 




6.1 


5.4 


16.0 


21.9 


T - 


ep^ 


6.7 


5.7 


15.6 


25.1 


T - 


/^</' 


3.7 


3.8 


14.2 


19.9 


T 


e(j) 


4.1 


4.5 


14.0 


22.0 


r - 


-> yua; 


7.0 


8.9 


25.7 


29.0 


r - 


-)■ euj 


8.6 


9.7 


21.1 


37.1 


r — 


' pK*^ 


4.9 


5.2 


15.8 


21.2 


T — 


> eK*^ 


5.7 


6.7 


15.6 


25.1 


T — 


' pK*^ 


4.9 


5.2 


15.8 


21.3 


T — 


> ei?*o 


5.2 


5.7 


15.6 


24.6 



of a fake muon and two pions. For the r — jj^uj mode, there is a small background in the 
signal region from generic t~^t~ events in which a fake muon is combined with two pions 
from r~ — tt^tt'^it^i't- and a fake tt" (from ISR and beam background). A large background 
is from — h^oouT in the AE < GeV and M^yo < region. For the r — )■ /i0 mode, 
the dominant background is due to qq and t~^t~ events in which a pion is misidentified as 
a kaon. For the r — )■ ^K*^ and /x^*'' modes, the dominant background comes from generic 
r+r~ events with a fake muon, fake kaon and real pion from — > 7T~it~^7t~i't-. If one of the 
pions is misidentified as a kaon, then the reconstructed mass from generic r+r" background 
could be greater than the r lepton mass, since the kaon mass is greater than that of the 
pion. 

For the fip^ and fiu modes, where significant background remains after applying the 
above-mentioned selection criteria, we extrapolate to the signal region by fitting to observed 
data in the ±5(Ta£; M^yo data sideband using the sum of an exponential and a first-order 
polynomial function for generic t~^t~ and continuum, respectively. For the r — )■ eV^, p,(j), 
fiK*^ and ^K*^ modes, the background after all the event selection requirements have been 
applied is small; extrapolation to the signal region assumes that the background distribution 
is fiat along the M^yo axis within iSaAs- We estimate the expected number of the back- 
ground events in the signal region for each mode using the number of data events observed 
in the sideband region inside the horizontal lines but outside the signal region. The signal 
efficiency and background level for each mode are given in Table lllli After estimating the 
background, we open the blinded regions. We observe one candidate event for each of the 
r — )■ piK*^ and piK*^ modes, and no candidate events for other modes. In each case the 
number of events observed in the signal region is consistent with the expected number of 
background events. 

The dominant systematic uncertainties in this analysis relate to tracking efficiencies and 
particle identification. The uncertainty due to the track finding efficiency is estimated to 
be 1.0% per charged track, therefore, the total uncertainty due to the charged track finding 
is 4.0%. The uncertainty due to particle identification is 2.2%, 1.9%, 1.3% and 1.8% per 
electron, muon, pion and kaon, respectively. The uncertainties due to the V'^ branching 
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FIG. 4: Scatter-plots in the ±20a area of the M^yo - AE plane: (a), (b), (c), (d) and (e) show 
the np^, fi(j), nio, fiK*'^ and iJ,K*^ modes, respectively. Experimental data events are shown as 
solid circles. The filled boxes show the MC signal distributions with arbitrary normalization. 
The elliptical signal (blind) regions shown by solid curves are used for evaluating the signal yield. 
The region between the horizontal solid lines (excluding the signal region) is used to estimate the 
expected background in the elliptical region. 



fractions are estimated to be 1.2% and 0.7% for the (p and u, respectively. The uncertainty 
due to the integrated luminosity is estimated to be 1.4%. The uncertainties due to the 
trigger efficiency and MC statistics are negligible compared with the other uncertainties. 
All these uncertainties are added in quadrature, and the total systematic uncertainties for 
the different modes range from 5.3% to 6.2%, as shown in Table UTTl 



UPPER LIMITS ON THE BRANCHING FRACTIONS 

Since no statistically significant excess of data over the expected background is observed 
in the signal region, we set upper limits on the branching fractions based on a frequentist 
approach [2l[. We calculate the 90% CL. upper limit on the number of signal events 
including systematic uncertainty (sgo) from the number of expected background events, the 
number of observed events and the systematic uncertainty, using the POLE program without 
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FIG. 5: Scatter-plots in the ±20a area of the Mgyo - AE plane: (a), (b), (c), (d) and (e) show 
the ep^, ecp, euj, eK*^ and eK*^ modes, respectively. Experimental data events are shown as 
solid circles. The filled boxes show the MC signal distributions with arbitrary normalization. 
The elliptical signal (blind) regions shown by solid curves are used for evaluating the signal yield. 
The region between the horizontal solid lines (excluding the signal region) is used to estimate the 
expected background in the elliptical region. 



conditioning 22|]. The upper limit on the branching fraction [B) is then given by 



where A',-,- is the number of r~'^r~pairs, and e is the signal efficiency including the branching 
fraction of . The value A^^-.^ = 782 x 10^ is obtained from the inte grat ed luminosity times 
the cross section of r-pair production, which is calculated in KKMC [23| to be a-^r = 0.919 ± 
0.003 nb and = 0.875±0.003 nb for 782 fb'^ at T(4^) and 72 fb'^ at T(5^), respectively. 
Table UTTl summarizes information about upper limits for all modes. We obtain the following 
ranges for 90% CL. upper limits on the branching fractions: B{t — eV^) < (1.8—4.8) x 10~^ 
and B{t — > P'V^) < (1-2 — 8.4) x 10"*^. These results improve upon our previously published 
upper limits [lO| by factors of up to 5.7. They are also more restrictive than the recent 
results from BaBar (ill . I2I . 
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TABLE III: The signal efficiency (e), the number of expected background events (A^bg) estimated 
from the sideband data, total systematic uncertainty (usyst)) the number of observed events in 
the signal region (Aobs)) 90% C.L. upper limit on the number of signal events including system- 
atic uncertainties (sgo), 90% C.L. upper limit on the observed branching fraction (;Bobs.) for each 
individual mode. 



Mode 


e (%) 


Nbg 


O-syst (%) Aobs 




Bohs (xlO-8) 


r~ 




7.09 


1 AO 1 o o r 

1.48 ± 0.35 


5.3 





1.34 


1.2 


r~ 




7.58 


0.29 ±0.15 


5.4 





2.17 


1.8 




- fi-4> 


3.21 


0.06 ±0.06 


5.8 


1 


4.24 


8.4 


T' 


" -> e~(j) 


4.18 


0.47 ±0.19 


5.9 





2.02 


3.1 


t" 




2.38 


0.72 ±0.18 


6.1 





1.76 


4.7 


t' 


~ — )• e^uj 


2.92 


0.30 ±0.14 


6.2 





2.19 


4.8 


r~ 


p"K*^ 


3.39 


0.53 ±0.20 


5.5 


1 


3.81 


7.2 


r~ 


e-K*° 


4.37 


0.29 ±0.14 


5.6 





2.17 


3.2 


T~ 


p-R*^ 


3.60 


0.45 ±0.17 


5.5 


1 


3.90 


7.0 




e~K*° 


4.41 


0.08 ±0.08 


5.6 





2.34 


3.4 



SUMMARY 

We have searched for lepton-fiavor-violating r decays into a lepton and a vector meson 
using 854 fb~^ of data collected with the Belle detector at the KEKB asymmetric-energy 
e~^e~ collider. Since no evidence for a signal is found, we set the following 90% C.L. upper 
limits on the branching fractions: B{t — > eV^) < (1.8 — 4.8) x 10~^ and B{t — )■ /iV^°) < 
(1.2 — 8.4) X 10~^. These results improve upon our previously published upper limits by 
factors of up to 5.7. This improvement results both from using a larger data sample and 
from an improved rejection of specific backgrounds, such as di-baryon production in the 
continuum for the r — >■ ^V^ modes, and r~ — h^n^Vj. decays with a photon conversion for 
the r — >■ eV^ modes. The more stringent upper limits reported here can be used to constrain 
the parameter space in various models of new physics. 

Acknowledgments 

We thank the KEKB group for the excellent operation of the accelerator, the KEK cryo- 
genics group for the efficient operation of the solenoid, and the KEK computer group and the 
National Institute of Informatics for valuable computing and SINET3 network support. We 
acknowledge support from the Ministry of Education, Culture, Sports, Science, and Tech- 
nology (MEXT) of Japan, the Japan Society for the Promotion of Science (JSPS), and the 
Tau-Lepton Physics Research Center of Nagoya University; the Australian Research Council 
and the Australian Department of Industry, Innovation, Science and Research; the National 
Natural Science Foundation of China under contract No. 10575109, 10775142, 10875115 
and 10825524; the Ministry of Education, Youth and Sports of the Czech Republic under 
contract No. LA10033 and MSM0021620859; the Department of Science and Technology 
of India; the BK21 and WCU program of the Ministry Education Science and Technology, 



12 



National Research Foundation of Korea, and NSDC of the Korea Institute of Science and 
Technology Information; the Polish Ministry of Science and Higher Education; the Ministry 
of Education and Science of the Russian Federation and the Russian Federal Agency for 
Atomic Energy; the Slovenian Research Agency; the Swiss National Science Foundation; 
the National Science Council and the Ministry of Education of Taiwan; and the U.S. De- 
partment of Energy. This work is supported by a Grant-in-Aid from MEXT for Science 
Research in a Priority Area ("New Development of Flavor Physics"), and from JSPS for 
Creative Scientific Research ( "Evolution of Tau-lepton Physics" ) . 



[1] J. E. Kim, P. Ko and D. G. Lee, Phys. Rev. D 56, 100 (1997). 

[2] A. Ilakovac, Phys. Rev. D 62, 036010 (2000). 

[3] D. Black et al, Phys. Rev. D 66, 053002 (2002). 

[4] C.-H. Chen and C.-Q. Ceng, Phys. Rev. D 74, 035010 (2006). 

[5] E. Arganda et al, JHEP 0806, 079 (2008). 

[6] R. Benbrik and C. H. Chen, Phys. Lett. B 672, 172 (2009). 

[7] Z. H. Li, Y. Li and H. X. Xu, Phys. Lett. B 677, 150 (2009). 

[8] A. Abashian et al. (Belle Collaboration), Nucl. Instr. and Meth. A 479, 117 (2002). 

[9] S. Kurokawa and E. Kikutani, Nucl. Instr. and Meth. A 499, 1 (2003), and other papers 
included in this Volume. 

[10] Y. Nishio et al. (Belle Collaboration), Phys. Lett. B 664, 35 (2008). 

[11] B. Aubert et al. (BaBar Collaboration), Phys. Rev. Lett. 103, 021801 (2009). 

[12] B. Aubert et al. (BaBar Collaboration), Phys. Rev. Lett. 100, 071802 (2008). 

[13] K. Hanagaki et al, Nucl. Instr. and Meth. A 485, 490 (2002). 

[14] A. Abashian et al, Nucl. Instr. and Meth. A 491, 69 (2002). 

[15] S. Jadach et al., Comp. Phys. Commun. 130, 260 (2000). 

[16] D. J. Lange, Nucl. Instr. and Meth. A 462, 152 (2001). 

[17] S. Jadach et al, Comp. Phys. Commun. 70, 305 (1992). 

[18] F. A. Berends et al, Comp. Phys. Commun. 40, 285 (1986). 

[19] S. Brandt et al, Phys. Lett. 12, 57 (1964); E. Farhi, Phys. Rev. Lett. 39, 1587 (1977). 

[20] K. Hayasaka et al (Belle Collaboration), Phys. Lett. B 613, 20 (2005). 

[21] G. J. Feldman and R. D. Cousins, Phys. Rev. D 57, 3873 (1998). 

[22] See http://www3.tsl.uu.se/~conrad/pole.html; J. Conrad et al, Phys. Rev. D 67, 012002 
(2003). 

[23] S. Banerjee et al, Phys. Rev. D 77, 054012 (2008). 



13 



